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In living cells, proteins combine 3D bulk diffusion and ID sliding along the DNA to reach a 
target faster. This process is known as facilitated diffusion, and we investigate its dynamics in 
the physiologically relevant case of confined DNA. The confining geometry and DNA elasticity are 
key parameters: we find that facilitated diffusion is most efficient inside an isotropic volume, and 
on a flexible polymer. By considering the typical copy numbers of proteins in vivo, we show that 
the speedup due to sliding becomes insensitive to fine tuning of parameters, rendering facilitated 
diffusion a robust mechanism to speed up intracellular diffusion-limited reactions. The parameter 
range we focus on is relevant for in vitro systems and for facilitated diffusion on yeast chromatin. 
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I. INTRODUCTION 

The binding of proteins to DNA or chromatin is at 
the basis of gene regulation, as for instance transcription 
factors stick to a gene along the DNA to promote or 
inhibit the subsequent binding of RNA polymerases, and 
hence to down or upregulate its transcription . 

In 1970 Riggs found that the association rate of the 
/ac-repressor in E.coU is 10^^M~^5~^, two orders of mag- 
nitude larger than it would be if the reaction time were 
dominated by the 3-dimensional (bulk) diffusion of the 
proteins through cytoplasm [4]. How can proteins move 
around so quickly in the cell? An intriguing mechanism, 
alternative to bulk diffusion, was proposed by Berg and 
von Hippel in [5 . They conjectured that proteins may 
bind non-specifically to DNA and diffuse along it - the 
linear portion Is of DNA which is explored before the pro- 
tein detaches from it is known as the sliding length. DNA 
looping may also allow hops and jumps of the protein 
between different sites on the DNA [6l [7] . The expres- 
sion facilitated diffusion indicates all the different mech- 
anisms through which proteins exploit the presence of 
DNA in the cell in order to reduce the mean search time 
(hereon MST), namely the average time that one pro- 
tein needs to reach its target. By considering a protein 
performing a very large number of slides, attachments 
and detachments, it is possible to theoretically estimate 
the value of the MST as a function of the sliding length 
Is (see e.g. [2 [3| 18j). The theory leads to the intrigu- 
ing result that there should be an optimal sliding length 
for which the mean search time is minimum, and sig- 
nificantly smaller than in the bulk diffusion case, which 
provides an explanation for Riggs' experimental results. 

There are to date very few direct dynamical numer- 
ical simulations of facilitated diffusion (for an up-to- 
date theoretical review see [2 ). Ref. [9 simulated de- 
tailed protein-DNA interactions with molecular dynam- 
ics, focussing on a small, unconfined, DNA loop, while 
Refs. pi)-14 focused on few selected isotropic geome- 
tries and on a single protein diffusing inside the cell. Our 
aim here is to extend the study of facilitated diffusion 



in a number of ways, with a view to bridge the gap to- 
wards more realistic conditions. We start from the ob- 
servation that DNA is strongly confined in vivo [15 , and 
consider both isotropic and anisotropic containers, which 
are relevant for eukaryotic nuclei and bacterial cells re- 
spectively. Secondly, to assess the role of DNA confor- 
mations, we simulate different persistence lengths: this 
is important in eukaryotes as the persistence length of 
chromatin varies in the 40-200 nm range and is thought 
to be activity-dependent [16]. Finally, we extend our 
investigation to the case in which (9(10) identical pro- 
teins simultaneously search for their target - this is a 
more realistic scenario given e.g. the typical copy num- 
bers of transcription factors inside a bacterium such as 
E. coli [171 [18]. As we shall show, each of these realistic 
add-ons leads to important qualitative differences in the 
physics of facilitated diffusion. 



Specifically, we find that the confining geometry and 
the DNA elasticity affect the mean search time in a ma- 
jor way: our simulations show that anisotropy slows down 
facilitated diffusion, while flexible conformations lead to 
a larger speed-up around the optimal sliding length with 
respect to semiflexible ones. These differences are much 
more pronounced for short chains. Our main result is 
that considering a cell with multiple copies of a protein 
looking for their target, as is the case in nature, leads 
to a very different view on facilitated diffusion with re- 
spect to the standard one, based on single protein models. 
Instead of displaying a well-defined minimum as a func- 
tion of the sliding length, the MST for multiple proteins 
(i.e. the time one protein first hits the target) is remark- 
ably flat, so that facilitated diffusion offers a robust way 
to speed up the search, not relying on any optimisation 
linked to fine tuning of the parameters controlling the 
intracellular protein dynamics. This trend holds for the 
chain lengths considered here, which are more than one 
order of magnitude larger than the sliding length: how- 
ever, very long chains, which we have not simulated, may 
again display a different behaviour. 
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II. METHOD 

Our approach is based on direct Monte-Carlo simu- 
lations of the whole facilitated diffusion process. For 
the DNA, we use a coarse-grained bead-spring model 
where beads are joined by FENE springs and interact 
via a shifted and cut-off Lennard- Jones (LJ) potential. 
This provides a hard-core for the beads, hence captur- 
ing the polymer self- avoidance. The cut-off is at 26 a, 
while the energy scale associated with the LJ potential is 
^ = ^ksT. Within this setting the natural length scale in 
our simulations is provided by the size of the chain beads, 
a. We chose to simulate confining geometries with sizes 
of ~ lOcr - 30<j as detailed below, and chains with con- 
tour length Lc ranging from 60a to 200cr. Such relatively 
short polymers allow us to more stringently test the the- 
ories which typically rely on the existence of many un- 
binding and binding events and hence implicitly assume 
long polymers. Moreover conditions similar to our sim- 
ulations may be designed experimentally by using naked 
DNA {(J = 2.5 nm) of length 150-500 nm confined in vesi- 
cles of size 50 nm - a set-up which should be realisable 
with single molecule experiments in vitro. An intriguing 
alternative mapping is to yeast chromatin in vivo^ where 
the appropriate coarse graining entails a value of a ^ 30 
nm (containing ~ 3 kilo-base pairs of DNA), and typical 
chromosome length is ~ 10^ base pairs (or Lc ~ 300cr). 

The elasticity of the DNA is modelled by a standard 
Kratky-Porod potential, Ub — —K^j^.ti • t^+i, where 
ti = r^+i — r^, fi being the position of the i-th monomer, 
and where the persistence length is determined in terms 
of the bending rigidity, via Ip = Ki^cj /ksT . We have 
considered the case of semifiexible {Ip = 20<j) and flexi- 
ble (i^5 = 0) polymers, relevant to DNA and yeast active 
chromatin respectively [16 . Equilibrium DNA configura- 
tions are grown and relaxed inside cylindrical or spherical 
containers. For the confining geometry, we used cylinders 
having a total length of 30<j and radii Rcyi going from 4(j 
to lOcr to provide different aspect ratios for the container. 
Sphere radii Rs range from 7a to 13cr giving isotropic con- 
fining regions with the same volumes of the cylindrical 
counterparts. Rather than explicitly modelling polymer 
dynamics, the equilibrium DNA conformations are kept 
frozen during single diffusion processes. This quenched 
approximation for the DNA is in line with most previ- 
ous theoretical and numerical work [8], [10] and can be 
justified by noting that the typical time needed by the 
polymer to change from one conformation to a completely 
uncorrelated one (Rouse time tr [19 ) is ~ L^, and for 
the cases we consider is larger than the timescale associ- 
ated with protein diffusion, ^ R'^/Ds. 

Proteins are modelled as pointlike and non-mutually- 
interacting particles, whose initial positions are chosen 
to be uniformly distributed at the walls of the confining 
geometry. In the bulk, proteins diffuse with a coefficient 
Ds = 0.021cr^/ ns. We use reflecting boundary condi- 
tions at the surface of the sphere or cylinder. Interaction 
with DNA is modelled as follows: whenever the new pro- 




FIG. 1: (Color online) Some simulation snapshots. A semi- 
flexible polymer is confined inside a cylinder (a) or a sphere 
(b) , while some proteins (red dots) diffuse in 3D or slide along 
it. 

posed position for the protein falls inside one of the DNA 
beads, this starts sliding with probability p. During each 
time step spent sliding, the protein may hop either to the 
right or to the left DNA bead, or detach from the genome 
with a probability 1 —p. Different choices of p correspond 
to different values of the sliding length /g, the key opti- 
misation parameter in the theory of facilitated diffusion. 
The duration of a sliding step, Ats = 48 ns, was chosen 
in order to give a one-dimensional diffusion coefficient 
Di equal to in agreement with previous theoretical 
estimates and simulations [TO] (even though most ex- 
perimental results suggest Di < [HI [20]). When 
the particle is at one of the two terminal beads of the 
chain it is always reflected, as in [10 . Some snapshots of 
our simulations are shown in Fig. la and lb respectively 
for cylindrical and spherical confinement for semifiexi- 
ble chains. The diffusion process for one protein stops as 
soon as it reaches the target, here modelled as an internal 
DNA bead. The interaction radius between the protein 
and the target, was equal to r5 = 0.35cr (which corre- 
sponds to about 1 nm for facilitated diffusion on naked 
DNA). For a fixed DNA equilibrium configuration the 
MST is estimated by averaging the time needed by one 
particle to reach the target over 100 initial conditions. 
A second average is then performed over ~ 50 different 
DNA configurations equilibrated in the confined region. 



III. RESULTS 

A. Comparison with theoretical results 

We first compare our results to existing theoretical pre- 
dictions by focusing on the cylindrical confining geome- 
try and on the facilitated diffusion of a single protein. 
Empty circles in Fig. 2 correspond to the simulated de- 
pendence of the MST on Is in the case of Ip = 20cr, and 
Lc = 200<j (error bars are smaller than symbols). These 
points are compared with theoretical curves (solid lines) 
predicted in [8 (tm) and in [10 {tk)- The latter gives a 
better estimate of the MST for small Ig , as well as of the 
optimal sliding length. For large /g, according to both 
theoretical models, the MST should increase much faster 
than in our data, which instead tend to quickly saturate 
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FIG. 2: (Color online) Dependence of the MST on Is is the 
case of a cylinder having radius Rcyi — 9cr and Lc — 200cr 
according to our data and to the theoretical predictions made 
in 8 (tm) and in [10 (tk). The case shown in the inset 
refers to short chains (Lc = 60cr) and to very anisotropic 
confinement (cylinder with radius Rcyi = 4cr and side length 
L = 30cr). 
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FIG. 4: (Color online) The normalised histogram of the num- 
ber of particles which have reached the target at a given time 
t, in the case of A^p = 10 (blue, top curve) and Np = 1 
(red, bottom curve). Both graphs refer to a semiflexible poly- 
mer with Lc — 60cr and to a cylindrical confinement with 
Rcyi = 4cr. 
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FIG. 3: (Color online) Comparison between the dependence 
on Is of the MST for one protein, scaled by 10 (mean field 
approximation) and for the first of ten particles (numerical 
data). Data refer to semiflexible DNA in a cylinder {Rcyi — 
4a and Lc = 200cr). 



around a constant value. This discrepancy is more dra- 
matic for shorter DNAs (see inset of Fig. 2, corresponding 
to a volume fraction of ~ 2% - a similar albeit smaller 
effect is observed for spherical confinement, discussed be- 
low) where it is more likely for the diffusing protein to 
reach the target in just a few slides. Such extreme events 
violate the hypothesis, used to estimate the MST the- 
oretically, that the length scale over which the sliding 
mechanism occurs is much smaller than Lr [2T1. 



B. Search time for ten proteins 

So far we have considered the facilitated diffusion of a 
single protein. However in a cell there are typically more 



copies of the same protein which simultaneously search 
for the same target - and a reaction is initiated once the 
first reaches it. 

It is then natural to look at the time at which the 
first of several proteins finds the target. In Fig. 3 we 
plot the MST of the first of ten non-interacting and si- 
multaneously diffusing particles - for comparison we also 
show the corresponding MST curve for a single protein, 
rescaled by a factor of 10 as expected from a naive mean- 
field argument (critically analysed in Ref. [17], see also 
below). The simulated data largely deviate from this 
simple estimate. 

The ultimate reason for this deviation lies in the 
markedly non-exponential nature of the search time 
distributions, which suggests the existence of multiple 
timescales in the dynamics. This can be appreciated by 
plotting the distribution of arrival times, (i.e. the time to 
first reach the target), P(t), shown in Fig. 4 for both the 
case of one single particle and for that of ten particles. 
From the linear-log plot in Fig. 5 it is apparent that nei- 
ther distribution is exponential - a better fit is provided 
by a streched exponential distribution (data not shown). 

We note that, as our ten proteins are non-interacting, 
the probability that one of the Np diffusing particles 
reaches the target at time Rn {t) = Jq P{t')dt' ^ is given 

by 



RN,{t) = l-{l-R{t)) 



(1) 



where R{t) is the corresponding single protein probability 
(see e.g. [I3]). While directly considering the trajectories 
of 10 proteins in our simulations, we have verified that 
this formula holds for our data. If P{t) were a simple ex- 
ponential, this formula shows that there would simply be 
a factor of Np between the MST of a single diffusing par- 
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FIG. 5: (Color online) These graphs show the same data as 
Fig. 4 (Lc — 60cr, semiflexible configurations confined in a 
cylinder with Rcyi — 4cr), but in a log-linear plot to stress 
the non exponential behavior of the distributions and they 
refer to the case of one particle (left) and of ten independent 
particles diffusing simultaneously (right). 



tide and the one where are simultaneously diffusing 
in the system. 

Importantly, we also find that the Ig dependence of 
the MST of the first of ten proteins does not display a 
clear minimum anymore. This suggests that, when more 
than one protein are searching for the target, there is 
effectively no optimal sliding length! In other words, 
while facilitated diffusion is still faster than bulk diffu- 
sion, there is no need for the protein to fine tune the 
sliding length on the DNA to optimize the search pro- 
cess; rather, the advantages due to facilitated diffusion 
are robust and largely independent of the protein-DNA 
affinity. While this holds up to the longer chain length 
simulated here, we cannot rule out that very long chains 
may display a different behaviour. 



C. Effect of confining geometry and DNA elasticity 

To further bridge the gap between simulations and bi- 
ologically relevant conditions, it is important to address 
the effect of the aspect ratio of the confining geometry 
and of DNA elasticity. As mentioned in the introduc- 
tion, there is a range of naturally occurring values for 
these parameters, and it is interesting to explore what 
impact tuning these has on the dynamics of facilitated 
diffusion. For instance, is there a difference in the ef- 
ficiency of facilitated diffusion within eukaryotic nuclei, 
which are roughly spherical, and inside elongated bac- 
terial cells? Do proteins search faster for a target inside 
euchromatin regions of the chromosomes, which are more 
flexible than inactive, heterochromatic ones [TF? To ad- 
dress these questions, we systematically analysed differ- 
ent confining conditions and polymer elasticities, and re- 
port below some typical results for the MST as a function 
of Is. 

We have considered ten different cylindrical geome- 
tries, of aspect ratio L/(2i?cyi) {L and Rcy\ are cylinder 
height and radius respectively, L is kept fixed) variable 
between 3/2 and 15/4 - we only report results on the 
MST for the most anisotropic case, as the difference with 
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FIG. 6: (Color onhne) Dependence of the MST on h for dif- 
ferent geometries and DNA elasticity. Data refer to a sin- 
gle protein, Rcyi = 4cr and correspondingly Rs = 7cr, when 
Lc = 60cr (a) and Lc = 200cr (b). (c) shows the MST for the 
first of ten particles for the same values of Rs and Rcyi , when 
Lc 200cr. 



respect to the case of spherical confinement is more read- 
ily appreciable there. As anticipated, we also compare 
the cases of flexible {Ip = 0) and semiflexible {Ip = 20 a) 
polymers. In the Appendix we provide more details on 
the conformations of the polymer on which proteins per- 
form their search for the target, and in particular discuss 
how they change with elasticity and confining geometry. 
Here we focus on the results for the MST. Each panel of 
Fig. [6] (from a to c) shows the Is dependence of the MST 
for four different cases: semiflexible and flexible DNA 
inside a cylinder, and semiflexible and flexible DNA in- 
side a sphere with the same volume. Note that to better 
highlight the effect of anisotropy on the MST, we focused 
on the most anisotropic cylindrical geometry, namely the 
one with aspect ratio L/2Rcyi = 15/4. 

For the smallest chains considered {Lc = 60<j, Fig. 6a) 
we observe major effects both of geometry and polymer 
elasticity, especially for small Is ^ 2a. In particular, in 
constrast to what happens when Ip = 20<j in a cylinder 
(Fig. 2 inset), in all other cases a well-defined optimal 
sliding length exists, and, more importantly, the MST is 
much smaller. This points out how both isotropy and 
flexibility can play an important role in enhancing the 
facilitated diffusion. Our results suggest that, at least 
for short chains, the search is faster within a sphere and 
on a flexible polymer. 

Fig. (6)3 shows the same plots as in Fig. [6^, but rela- 
tive to the case of Lc = 200cr. Here the effects of the 
confining geometry and DNA conformations are less sig- 
nificant. This is due to the fact that the volume fraction 
occupied by the DNA now increases, so that even the 
stiffest conformations develop loops and local contacts. 



5 



as they need to bend several times to fit in. Moreover, 
long polymers are more homogeneously distributed in the 
cell than short ones, which renders the local environment 
of a DNA segment less sensitive to the overall confining 
geometry, or polymer elasticity. Finally, Fig. ^ shows 
the MST for the first of ten proteins to reach the target, 
for Lc = 200(7. The data confirm the lack of an optimal 
sliding length. Interestingly, we also find that, as for the 
case of a single diffusing protein, the search is quicker 
within a sphere. The effect of fiexibility is still present 
but milder than for a single protein. 

An analysis of the statistics of 3D excursions in be- 
tween slides, given below, sheds more light on these re- 
sults. 



D. Statistics of jumps and hops 

Here we provide an analysis of how the occurence 
probability of different SD-diffusion mechanisms such as 
"jumps" and "hops" is affected by DNA elasticity and 
confinement geometry, with a view of explaining our re- 
sults in Fig. [6] on the infiuence of these factors on the 
MST. 

"Hops" and "jumps" are 3D-diffusion events that take 
place between two slides. In a hop the distance along 
the DNA between the take-off (detachment) and land- 
ing (re-attachment) point is "small", while in a jump it 
is "large". In our treatment we consider the distance 
"small" if it is less than 10 beads, and "large" otherwise. 

The key features of the statistics of jumps and hops 
are shown in Fig. jT] The absolute numbers of jumps and 
hops give the clearest signal. As these quantities depend 
on the sliding length, we plot their value at the minimum 
MST (as a function of Is) in Fig. [7^,b. For Lc = 60cr, 
both quantities are largest in the case of a semifiexible 
DNA in cylindrical confinement, which correlates well to 
a larger MST in this situation. Furthermore, still for 
Lc = 60<j, in the case of spherical confinement a smaller 
number of jumps is required on a fiexible DNA (case 4 in 
Fig. [7| than on a semifiexible one (case 2 in Fig. [t]), to 
reach the target, again in agreement with the MST trends 
in Fig. [6j Facilitated diffusion on a fiexible DNA inside a 
cylindrical container would require the smallest number 
of jumps, which does not match the fact that the MST in 
this case is not optimal. This is because in this situation 
the number and duration of hops also plays a role: these 
are relatively numerous, and, more importantly, each of 
these lasts about 50% longer, with respect to the average 
of the other cases (see Fig. [Tji). 

It is also interesting to monitor the ratio between the 
number of jumps and the number of hops, which instead 
does not depend on the sliding length. This ratio is once 
again affected by the polymer elasticity, more so in the 
case of a cylindrical confinement, where for a semifiexible 
DNA it is equal to 0.11, almost twice as large as for a 
fiexible DNA (0.06). In the case of spherical confinement 
the ratio between hops and jumps is again equal to 0.11 
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FIG. 7: (Color online) Histograms describing the average 
number (a, b) duration (c, d) and the average distance be- 
tween the takeoff and landing points ('length', e, f) for jumps 
and hops. Blue (left) columns refer to Lc = 60cr, while orange 
(right) columns refer to Lc — 200cr. Numbers on the horizon- 
tal axis refer to different geometries and elasticities: (1) and 
(2) are for semifiexible DNA confined in a cylinder and in 
a sphere respectively, (3) and (4) refer to flexible DNA in a 
cylinder and in a sphere respectively. All data refer to the 
case of Rcyi — 4cr (L/(2Rcyi) = 15/4) - and correspondently 
Rs = 7a. 



for semifiexible polymers, and rises to 0.09 for fiexible 
ones. These differences are related to the average back- 
bone distance between the takeoff and landing points, in 
both jumps and hops: as can be seen in Fig. [7^ and f, 
particles perform longer jumps along a semiffexible than 
along a flexible polymer. In contrast, hops are shorter 
when performed along semiflexible configurations. This 
is evidence of the fact that stiffness favours either very 
short or very long range (in terms of chemical distance) 
three-dimensional displacements, while ffexibility allows 
intermediate-range displacements, due to enhanced lo- 
cal bending. Short-range 3D displacements are advanta- 
geous for short DNA chains, as witnessed by our MST 
data in Fig. [6) 

When longer polymers are considered {Lc = 200cr), 
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though, the trends are less clearcut. On the one hand, for 
very low values of /g, the MST in a cylinder is again larger 
with respect to that in a sphere, and once more this is in 
line with the trends in the total number of jumps. On the 
other hand, now the semiflexible DNA inside cylindrical 
confinement no longer gives the slowest search, especially 
for larger values of Z^, where the MST is largest for a 
fiexible chain inside a cylinder. This is in part because 
the frequent intermediate-range displacements performed 
on fiexible polymers confined in cylinders make a less 
efiicient strategy in longer chains, as this wastes time in 
sampling nearby (in terms of chemical distance) DNA 
regions. At the same time, it is interesting to note that 
jumps always take (slightly) longer along fiexible than 
along semifiexible chains. We ascribe this to the fact 
that in fiexible polymers most of the nearby beads are 
also adjacent along the chain, and it is less likely for a 
protein to find a closeby bead which is very far along 
the DNA. This is a disadvantage especially in the case 
of long polymers, where jumps are expected to play an 
important role in favouring a better sampling of DNA. 
As a result the MST inside a cylinder is slightly longer 
in the fiexible case. 

We finally note that in all cases isotropic confinement 
reduces the MST, for both small and long polymers. This 
can be explained by noting that in a sphere any two 
points are on average closer together than inside a cylin- 
der. This is quantified via Fig. 11 (see Appendix), show- 
ing the mutual distance distribution function, n(r) [22 . 
This distribution is much more spread for cylinders than 
for spheres, which allows a particle in a sphere to reach 
its target after a considerably smaller number of 3D- 
diffusion events. 

To summarise, in this section we have analysed the 
statistics of jumps and hops with a view of understand- 
ing more deeply the results in Fig. |6] on the effect of 
geometry and elasticity on the MST trends. Firstly, we 
have seen that, especially for short polymers, the MST 
trends can be rationalised by looking at the total num- 
ber of jumps and hops needed to get to the target. This 
may be viewed as a consequence of the fact that the mu- 
tual distance distribution function goes to zero at a much 
shorter distance in an isotropic geometry. Secondly, we 
have shown that in a cylinder, 3D-displacements along a 
fiexible chain usually take longer than in the presence of 
semifiexible DNA. This increases the MST of long poly- 
mers, where the search requires several such jumps. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented detailed Monte Carlo 
direct dynamical simulations of the facilitated diffusion 
of one or several proteins inside a cylindrical or spheri- 
cal cell in presence of equilibrated DNA configurations. 
We have compared our simulation results with theoreti- 
cal predictions for the mean search time, the time needed 
for one protein to first reach its target, for instance a pro- 
moter along the genome. For short DNA chains we find 
a large discrepancy, which we attribute to the fact that 
proteins can reach the target with no or few sliding steps 
along the genome. This discrepancy attenuates as the 
contour length of the DNA increases. Our main result is 
that the average search time needed for the first of sev- 
eral copies of the same protein looking for their target si- 
multaneously shows a qualitatively different dependence 
on physical parameters: most notably, there is no longer 
an optimal sliding length leading to the fastest search. 
This is important biologically, as in bacterial or eukary- 
otic cells the typical copy numbers of e.g. transcription 
factors is low but larger than unity (typically '^10 [18 ). 
This fact may then be interpreted as a way to ensure that 
the search for the specific target of a protein is not cru- 
cially dependent on the precise value of its non-specific 
affinity with DNA. Interestingly, we have also found that 
polymer elasticity and the confining geometry may both 
affect the mean search time: facilitated diffusion is most 
effective inside spherical cells and in presence of flexible 
polymers. We flnd it intriguing to note that these advan- 
tages may be exploited for facilitated diffusion inside the 
spherical nuclei of eukaryotic cells such as yeast, where 
an enrichment in local gene activity is thought to corre- 
late with an enhanced chromatin flexibility [16]. Finally, 
we hope that our results may stimulate in vitro single 
molecule experiments on facilitated diffusion of proteins 
on DNA molecules of different contour lengths and under 
different conflning geometries. 
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FIG. 8: (Color online) Distribution of the distance of beads 
from the cylinder axis for semiflexible and flexible DNA (left) 
and of the distance from the cylinder and sphere centre for 
semiflexible DNA (right). Different colours (symbols) refer to 
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squares) and 200cr (green, triangles). Both graphs refer to 
Rcyi = 4(7, and correspondently to Rs = 7a. 
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Appendix A: Characterisation of DNA 
conformations 

In this Appendix we give more details on the DNA 
conformations, and on their dependence on the polymer 
flexibility and on the confining geometry. 

Fig. 8 (left) shows the distribution of the distance of 
DNA beads from the cylinder axis. Semiflexible polymers 
tend to concentrate more near the walls, since this allows 



FIG. 9: (Color online) Dependence of the initial diffusion 
time, spent on average by one particle in performing 3D- 
diffusion in the bulk before starting interacting with DNA, as 
a function of the volume fraction. Different colours (shades of 
gray) refer to different confining geometries and elasticities: 
semiflexible (black) and flexible (light gray) polymers inside 
a cylinder (filled symbols) and a sphere (open symbols). Val- 
ues for different DNA contour lengths are shown: Lc — 60a 
(triangles), Lc = 120cr (squares) and Lc — 200a (circles). 



them to minimise the bending. Flexible polymers instead 
are more homogeneously distributed inside the confining 
geometry. Fig. 8 (right) shows the distribution of the 
distance of the beads of a semiflexible DNA from the 
centre of the confining geometry (either a cylinder or a 
sphere), in the case of the most anisotropic cylindrical 
geometry {L/{2Rcyi) = 15/4) and of a sphere having 
the same volume {Rg = 7a). While in the case of a 
cylinder all distances between Rcyi = 4cr and L/2 = Iba 
are almost equally likely, in the sphere isotropy leads to 
a peak at distances slightly smaller than Rs. [Similar 
results are obtained when comparing the distribution of 
beads for semiflexible and flexible DNA in a sphere or for 
flexible DNA in different geometries.] 

The spatial distribution of the DNA correlates well 
with the "initial diffusion time", i.e. the time that a 
protein, starting from the surface of the cell (whether 
spherical or cylindrical in shape) takes before "touching" 
the DNA for the first time, independently of whether or 
not a slide is performed then. This initial diffusion time, 
computed for the different DNA conformations consid- 
ered, is shown as a function of the volume fraction in 
Fig. 9. This quantity is larger in the case of flexible 
polymers, and it is maximal for short flexible chains, in 
which case the initial diffusion time impacts considerably 
on the total MST. This trend, however interesting, is not 
impacting much on the MST, whose dependence on ge- 
ometry and flexibility is more crucially determined by the 
average number of jumps and hops (see main text). 

A more detailed characterisation of DNA elasticity 
is provided by the tangent-tangent (normalised) corre- 
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FIG. 10: (Color online) The graph shows the dependence of 
the persistence length on the volume fraction in the case of 
Lc — 200(7 and Rcyi — ^cr for semiflexible DNA confined in 
a cylinder. Different colours (shades of gray) refer to diffent 
confining geometries and elasticities: semiflexible (black) and 
flexible (light gray) polymers inside a cylinder (filled symbols) 
and a sphere (open symbols). Values for different DNA con- 
tour lengths are shown: Lc = 60cr (triangles), Lc = 120cr 
(squares) and Lc = 200cr (circles). 



lation function Ct{m) = {U • ti-^m) /\U\\U-\-m\^ where 
ti = r^+i — r^, and fi gives the position of the i-th 
bead. The characteristic decay length of Ct{m) (which is 
defined by assuming an exponential decay) defines the 
persistence length. In the unconfined case, for semi- 
flexible polymers this is Ip = 20<j as expected. When 
DNA is confined, though, it needs to bend several times 
to fit in the confining geometry, and Ct consequently 
decreases more rapidly in space, so that the measured 
effective persistence length, Ip^eff^ defined according to 

Ct{m) ~ e ^p^^ff is always smaller than Ip. Fig. 10 shows 
Ip^eff as a function of the volume fraction in the four dif- 
ferent cases of a semiflexible and of a flexible polymer 
in a cylinder and in a sphere. It is interesting to notice 
here that the isotropic confinement of the sphere deter- 
mines a smaller Ip^eff with respect to the one measured in 
a cylinder, for semiflexible configurations with the same 
contour length. This is due to the fact that the DNA, in 



order to fit into a sphere, needs to bend more than in a 
cylinder with the same volume. 

Finally, an important DNA configurational property 
affecting the MST is given by the "local environment" 
around a given bead of the chain that can be probed, 
for example, by the distribution probability of all mutual 
distances (between any two beads in the chain), which we 
call n(r). This is plotted in Fig. 11, only for beads whose 
distance along the chain is larger than 5<j (this threshold 
is chosen to filter out the information from consecutive 
beads along the chain). Fig. 11 shows that for Lc = 
60(7 the tvDical distance between two beads is larger in a 
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FIG. 11: (Color online) Dependence of the mutual distance 
distribution function, n(r) for different confinement geometry 
and polymer elasticity, see legend. Data correspond to Lc — 
60(7, to Rcyi = 4(7 and to Rs = 7(7. In the calculation only 
beads separated by 5(7 or more along the backbone of the 
polymer are included. 



cylinder, for semiflexible chains. This is because flexible 
chains are more compact on average so that points are on 
average closer together. Furthermore, in a sphere n(r) 
goes to zero at a smaller value of the interparticle distance 
since, despite having the same volume, cylinders are more 
elongated and therefore allow any two beads to be farther 
apart with respect to what happens in spheres, where 
the maximum distance between two points is set by the 
diameter. 



